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Abstract 

Optical lattices offer the possibility to investigate the superfluid properties of both Bose conden- 
sates and Fermionic superfluid gases. When a population imbalance is present in a Fermi mixture, 
this leads to frustration of the pairing, and the superfluid properties will be affected. In this con- 
tribution, the influence of imbalance on the flow of a Fermi superfluid through an optical lattice 
is investigated. The flow through the lattice is analysed by taking into account coupling between 
neighbouring layers of the optical lattice up to second order in the interlayer tunneling amplitude 
for single atoms. The critical velocity of flow through the lattice is shown to decrease monotonically 
to zero as the imbalance is increased to 100%. Closed-form analytical expressions are given for the 
tunneling contribution to the action and for the critical velocity as a function of the binding energy 
of pairs in the (quasi) two-dimensional Fermi superfluid and as a function of the imbalance. 
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I. INTRODUCTION 



In the past several years, ultracold Fermi gases have become a major focus of research 
on the interplay between Cooper pairing, strong interactions, and reduced dimensionality. 
The appeal of this system stems from the fact that both the interaction strength and the 
geometry of the confinement potential can be controlled very precisely. Feshbach resonances 
can be applied to tune the s-wave scattering length over a wide range of positive and negative 
values, which allows to investigate the crossover between a regime of Bose- Einstein condensed 
(BEC) molecules[l] and a Bardeen-Cooper-Schrieffer (BCS) superfluid [2)]. Optical lattices 
are used to reduce the dimensionality of the system to ID or 2D or to create a 3D optical 
lattice {3, 4 , sj] . Moreover, since the amount of atoms of each species in a Fermi mixture 
can also be controlled accurately, it has become possible to investigate the effect, on pair 



a. 



formation, of a population imbalance between the pairing partners [ 

In this contribution, we investigate the case of an imbalanced Fermi superfluid in a one- 
dimensional optical potential. In previous work, it was shown that the gap and number 



equations can be solved exactly for an imbalanced two-dimensional Fermi gas [13] - this 
corresponds to a single layer in the optical potential. Here, we will take into account the 
coupling between the different layers in the ID lattice, and derive expressions for the tun- 
neling amplitude and the critical velocity of the imbalanced Fermi superfluid through the 
optical lattice. Tunneling of a balanced Fermi superfluid through a ID optical lattice has 
been studied experimentally [5!], as has the imbalanced 3D Fermi superfluid jf], [7]. However, 
to the best of our knowledge, tunneling of a imbalanced superfluid has not yet been inves- 
tigated. 



II. OPTICAL LATTICE 



The system that we study in this contribution consists of a Fermi mixture with unequal 
amounts of two components, referred to as "spin up" and "spin down". This is loaded into 
a one-dimensional optical lattice (in the ^-direction). The optical lattice is created by two 
counterpropagating laser beams with wave length A and is modelled by a potential energy 
Vq sin 2 (2nz/X). Under the influence of this lattice, the gas forms a stack of typically a 
few hundred quasi-2D layers containing several thousands of atoms each. The confinement 
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influences the way that interactions are taken into account. The interaction between atoms 
within one quasi-2D layer can be modelled by a 2D-pseudopotential V = gS(r) where the 
strength g depends on the energy of the scattering atoms through 
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Here m is the mass of the atoms, and Ej, is the energy of the bound state that always exists 
in two dimensions 8[] , given by E b = (CHlul/tt) exp(v / 2vr^L/a s ) , with a s the (3D) s-wave 



scattering length of the fermionic atoms, ujl = ^/8tt 2 V / (m\ 2 ) and £l = yV ( muJ L) and 
C w 0.915 (cf. Ref. |9|). 

The goal of this contribution is to take into account the coupling between adjacent valleys 
of the optical potential. This coupling is due to tunne ling of individual atoms from layer to 
layer and is characterised by the tunneling amplitude [ll 
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III. PATH-INTEGRAL APPROACH 
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We follow the path-integral approach, as applied by Iskin and Sa de Melo [U[ for an 
imbalanced gas, and as applied in Ref. [12j to the optical lattice. The action functional S 
for the fermionic mixture in the optical lattice can be written as a path integral over the 
exponential of an action functional, consisting of the contributions Sj of the individual layers 
and the contributions Sjj+i of tunneling between adjacent layers. That is, we write for the 
action S = (Sj + Sj t j + %), with the single- layer action 

Si = EC H"» + k2 - ^) C + & { ^^%^V ( 3 ) 

k,a k 

and the tunneling action 

= E * + CC') • ( 4 ) 

k,cr 

Here, we use k = {k, u n } for the the 2D wave number k and the Matsubara frequency u n , and 
a for the spin, so that tjj^l and tp^l represent the Grassmann variables for the fermionic fields 
in layer j. The system parameters are the interaction strength g, the tunneling amplitude 
t, and the chemical potentials ^ that fix the amounts of spin-up and spin-down particles. 
We use units h = 2m = 1. 
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The partition sum is the functional integral of exp{— S} over the Grassmann variables. 
The interaction part is decoupled by introducing the Hubbard-Stratonovic fields Aj,Aj in 
each layer, after which the functional integral over the Grassmann variables can be taken. 
Writing those fields as a function of amplitude and phase, Aj = \Aj\ e j , we find the following 
effective action 

Seff = ~~ £ M " Tr {log [-Qji + T j>f ] } (5) 
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Here, the trace is taken over k, over the spin variables and over the layer index. The single 
layer inverse Green's function is 

-e- = ^xf- ! "" + k2 -"« |A)I ] ( 6 ) 

y \Aj\ -iu n -k 2 + fiu 

and the tunneling propagator is 

/ te ^-^)/2 o . 

%j> = i S f,j+l + S 3',i-l) x „ V9 ( 7 ) 

\ _t e »(*i-«j')/2 

Setting t = we obtain the result for the imbalanced 2D Fermi superfluid, discussed in 



Ref. [I3j|. We write this result as S e ff(t = 0) = £\ fiflj where (3 = l/k B T is the inverse 
temperature, and Qj is the BCS free energy for layer j. If, on the other hand, we keep t > 
but set fi^j = ui j w e obtain the results for the balanced Fermi gas in an optical potential, 
derived in Ref. [12]. Here, we keep t > and investigate fi^j ^ fi^. We now treat the 
additional tunneling terms as a perturbation, and expand up to order t 2 . Diagrammatically, 
the process we include in the self-energy of layer j consists of a particle tunneling back and 
forth over layer j + 1 or j — 1. The result is 



S eff ~ £ 



2t 2 cos(6> i+ i - 6j) \Aj\ \A j+1 \ 



(8) 



[(Q - iu n ) 2 - (4 J) ) 2 ][(0 + i - iunf - (4 J+1) ) 2 ] 

where J^ fc represents both the sum over Matsubara frequencies and the integral over k. Here, 
A*j = (a*t»j +A t ij)/2 is the average chemical potential and Q = (fi^j — //y)/2 is the difference 
in chemical potentials, and 



Ejp = y/(k*- N )' + \A j \ 2 , (9) 
^ = k 2 - N . (10) 



4 



Note that up to this order of perturbation, the normal particles do not contribute to tun- 
neling: when we interleave layers of superfluid with layers of normal gas (for example the 
excess spin population), the product \Aj\ |Aj +1 | vanishes. When we assume that no such 
interleaving takes place and that moreover |A 3 -| ~ |A J+1 |, the integration over k and the 
Matsubara summation can be performed analytically and results in 

S eff = Pi + T j+1J cos(% +1 - 9,)] , (11) 
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with 



7, +1J = il 1+ ^=- 2 Vi^,. (12) 



where for Q < |A,-| the last term is not present; this case corresponds to the known result 



for the balanced gas 



12|. 



IV. JOSEPHSON REGIME AND CRITICAL VELOCITY 

We will investigate the specific case where both the 2D density + n± and the imbalance 
do not change significantly over the layers. In that case expression (fT2"l) . is also 

not changing significantly from layer to layer, and we can drop the index j for |A| ,/x and 
(. The superfluid motion over the layers is then due to the phase differences over the 
layers. Using the results from [13J, we can re-express as a function of the imbalance 

5n/n = (rif — n^)/{n^ + n±), and of the binding energy per particle Ej,/(2Ep), in units of 
the Fermi energy Ep = (hkp) 2 /(2m) with kp = 2n(n^ + nj_). 

Three regimes can be identified. The first regime occurs for large imbalance, 5n/n > 
Eb/(2Ep). Then no superfluidity occurs - the Fermi mixture is too imbalanced to support 
pairing. The second regime is the 'weak pairing' regime, characterised by 5n/n < E^j (2Ep). 
In the three-dimensional case, one would refer to this as the BCS regime. Then we have, 
from Ref. [3], the following analytical results for the 2D case: 

|A| 2 = 2E b \l - y/2/E b (5n/n)\ (13) 
li = l- (E b /2) [l - y^2/E~ b (5n/n)j (14) 
C = ^/|A| 2 + (5n/n) 2 (15) 
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FIG. 1: Imbalance between the spin components reduces the critical velocity of a Fermi superfluid 
in an optical lattice. Here, the critical velocity relative to the critical velocity for a balanced gas 
is shown as a function of the pair binding energy E b and the imbalance Sn. The dashed line at 
E b = 2Ep separates the regime of weak pairing from that of strong pairing. For E b jEp < 25n/n 
(white region) superfluidity is suppressed. 



This results in 
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The third regime is that of strong pairing, characterised by Eb/(2Ep) > 1. In three 
dimensions this regime might be compared to the BEC regime. For the strong pairing 
regime, we find 13] 



|A| 2 = 2E b (l -dn/n) 

fi=l- (E b /2) (1 - 5n/n) 

C=J\A\ 2 + [2{5n/n)- h if 



(17) 
(18) 

(19) 
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which results in 




8tt 2E F + E b - y/2E^E b f 



t 2 AE F 



+2 



(2f-l)+f 



(20) 



^(2^-l) 2 +(2^ + l)^(l-^)+f (l-^) 2 



The presence of imbalance reduces from its result without imbalance. Since the 

critical velocity can be written as v c = XTj + i j/(hE F ) (see Ref. [yj]), we find that that 
the imbalance also reduces the critical velocity for flow through the lattice. In Fig. 1, 
we show this reduction v c (5n)/v c (5n = 0), as a function of E b and 5n/n. Note that for 
Eb/Ep < 2Sn/n superfluidity is not present. The two regimes of weak pairing and strong 
pairing can be distinguished by a kink in the contour lines. 

V. CONCLUSIONS 

The theory of flow through the lattice set up Ref. [3| is largely unmodified by the 
presence of imbalance. The main effect is the reduction of the tunneling coefficient 
expression (fl2l) . In this work, we kept 5n constant and assume that |A| , /i, ( vary slowly 
from one lattice site to another. Under these assumptions, we have derived closed analyt- 
ical formulae for the tunneling contributions (and thus the critical velocity) for flow of an 
imbalanced superfluid through an optical lattice. 

The assumption that | A| , fi,( vary slowly is good unless we have interleaving of normal 
gas and superfluid gas layers. If this is the case tunneling is suppressed because we take 
the optical potential to be deep enough such that the normal gas is pinned and flow is only 
due to phase coherence of the pair condensate. The assumption of constant 5n is related; it 
also relies on the fact that there is no phase separation. This will be much more difficult to 
satisfy in practice, as coexistence may only be achievable at finite temperatures. In practice, 
at low temperature, there will be phase separation in the layer and the tunneling will only 
take place in the region of superfluid phase, and will be suppressed in the normal layer 
around it. This dynamical interplay between pinned normal state, and a phase separated 
superfluid are outside the scope of this paper, which relies on the possibility to create (albeit 
dynamically) an imbalanced superfluid. 
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